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Abstract. Dirac electrons with a zero-gap state (ZGS) in organic conductor a-(BEDT- 
TTF) 2 l 3 result from a fine tuning of the seven nearest neighbors transfer integrals 
(ui, a 2 , 03 , 6 i, b 2 , fes, 64 ) between the four molecules of the unit cell. In this work we show that 
for given moduli |ai|, ...| 64 |, the possibility of having Dirac electron with a ZGS at 3/4 (or 1/4) 
filling strongly depends on the specific configurations of signs of the seven transfer integral. More 
precisely it is possible to classify the sign configurations into essentially four classes determined 
by Xa = sign(a 2 a 3 ) and Xb = sign(foi 62 fe 364 ). Using extended numerics, we show that for both 
weak and large inhomogeneity in the moduli, the class {Xa,Xb) = (~i~) is the most favorable 
to find Dirac electrons with ZGS at 3/4 (or 1/4) filling. For the class {Xa,Xb) = (+ 1 +) in the 
opposite case, we never found any ZGS at either 1/4 or 3/4 filling. The last two classes given by 
{Xa,Xb) = (+5 ~) and (xa, Xb) = (~i +) corresponding to an intermediate situation; they allow 
for ZGS at 3/4 (resp. 1/4) filling but are much less favorable than class (xa, Xb) = (~) “)• As a 
matter of fact, all previous numerical studies of Dirac electrons and ZGS in a-(BEDT-TTF) 2 l 3 
correspond to class {Xa,Xb) = (“!+)■ 


1. Introduction 

It is now well established that organic conductors q:-(BEDT-TTF) 2I3 [I], exhibit a zero gap- 
state (ZGS) characteristic of the presence of massless Dirac electrons close to the Fermi energy 
[ 2 ]. However, compared to graphene, the Dirac cones are strongly tilted resulting in a strong 
anisotropy of the Dirac electron velocity components. Furthermore, the fact that a-(BFDT- 
TTF)2l3 has a unit cell composed of four inequivalent molecules leads to some unique physical 
properties associated to this ZGS. In particular, as revealed by recent NMR measurements, the 
anisotropy of the Dirac particles appears strongly related to their local spectral weight on the 
four inequivalent molecules per unit cell; such that the smallest (largest) velocity Dirac particles 
have their largest weight on the electron rich (poor) molecule site O 01 [ 5 ] . 

From a theoretical point of view, the electronic properties of a-(BEDT-TTF)2l3 are usually 
described by a tight-binding model with seven nearest-neighbor transfer integrals oi,-- - ,64 
between the four molecules A, A’, B, and G of the unit cell. As illustrated on hgure Fig. this 
model presents an inversion symmetry with inversion centers located on B,G and the middle 
point between A and A’.|T]. Thank’s to this inversion symmetry, recent studies have used Fu- 
Kane topological argument in order to determine quantitatively the existence domain of Dirac 



Figure 1. (Color online) Structure of a-(BEDT-TTF) 2 l 3 , with four molecules A, A’, B, and C 
per unit cell (dotted square). Inversion centers are located on sites B, C, and the middle point 
between A and A’.[T]. The seven nearest-neighbor transfer integrals are shown by oi,--- , 64 . 
The chirality of four distinct triangular plaquettes is shown by xi, X 2 , XS; and X 4 . 


electrons from the knowledge of the energy and inversion parity eigenvalues at the four time 
reversal invariant momenta (TRIM) 13 IB]. Despite their usefulness, these studies do not 
allow to distinguish ZGS systems in which Dirac electrons are at the Fermi energy (with a 
vanishing density of states at the Fermi energy) from metallic systems in which there is an 
indirect overlapp between valence and conduction band such that the Dirac electrons are below 
(or above) the Fermi energy [9l [10] . 

The purpose of the present work is to clarify the relative role played by the the inhomogeneity 
in the moduli and in the signs of the seven transfer integrals on the appearance and stability of 
Dirac electrons with ZGS. In fact as explained in more details below, for given moduli |ai|, ...I 64 I, 
the possibility of having Dirac electron with a ZGS at 3/4 (or 1/4) filling strongly depends on 
the specific configurations of signs of the seven transfer integral. More precisely it is possible 
to classify the sign configurations into essentially four classes determined by Xa = sign(a 2 a 3 ) 
and Xb = sign( 6 i 62 & 3 ft 4 )- Using extensive numerics, we show that for weak, large and no 
inhomogeneity in the moduli, the class (xa, Xb) = (—, —) is most favorable to find Dirac electrons 
with ZGS at 3/4 (or 1/4) filling. The classes {xa,Xb) = (+)~) and {xa,Xb) = (~ 5 +) constitute 
equivalent but less favorable situation for the occurrence of Dirac electrons with ZGS at 3/4 (or 
1/4) filling. Finally, the last class is (Xa,Xb) = (+;+) and for the range of inhomogeneity that 
was explored no ZGS was obtained for either 1/4 or 3/4 filling. 

The paper is organized as follows. In §2 we present the model and the main quantities of 
interest. In §3 we first introduce the chirality patterns and their properties and then using 
properties established in Appendix A we explained the relation between the chirality patterns 
and the energy band structure and local density of states. In §4 we present our extensive 
numerical studies of the existence of ZGS at 3/4 filling for the three pertinent classes of chirality 
patterns. In §5 we summarized our main results. 

2. Tight-binding model 

The starting point is the 4x4 Bloch Hamiltonian matrix H(k) corresponding to the real space 
tight-binding model shown in Fig. 
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with a = as + 026*^^', 6 = 63 + 626®^^, c = e^^y{b4, + and / = oi(l + (the lattice 

constant is taken as unity). The Boch basis \jk) [j = A,A',B,C) was chosen such that the 
Hamiltonian matrix is Bloch periodic H{k + G) = H{k) with G a reciprocal lattice vector. 
In this Bloch basis, the energy eigenstate \En{k)) of the band is represented by a four 
components vector Un{kY' = ,u^,u^) such that H{k)un{k) = En{k)un{k) where by 

convention Ei{k) < E2{k) < E^lk) < E4^{k). The local density of states (LDOS) on each site is 
then obtained as (j = A, A', B, C) 


/ dk 

— Y^\ui{k)\H{E-En{k)). ( 2 ) 
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For 3/4 filling, the existence of a pair of Dirac points at Ekjy implies E^^kj:)) = E3{kD) = E^^ 

with ki) 7^ G/ 2 ‘, the existence of a ZGS further necessitates that E^^k) > Ej.^ > E^^k) for 

any k. In that situation the LDOS verify Dj{E) ~ \^-^d\ near Ed where vj is the site 

D 

dependent angular averaged Fermi velocity. 


3. Classification in terms of chirality patterns 

For fixed moduli (|ai|, • • • the energy bands En{k) and the eigenfunction components Un(k) 
have the freedom of the sign of the seven hopping par ameters, as a result one may expected 
to obtain 2 ^ distinct LDOS Dj{E). As explained in Appendix A, the number of distinct 


LDOS Dj{E) is entirely determined by the number of patterns of the four chiralities y* = ± 
(z = 1 , 2 , 3 , 4 ) [m 


Xi = sign(a 26 i 64 ), X2 = sign(a 362 & 3 ), Xs = sign(ai6i62), X 4 = sign(ai 6364 ), ( 3 ) 

where each chirality Xi encodes the sign of the product of the transfer integrals for each of the 
four distinct triangular plaquette that compose the unit cell. The are in total 2 “^ = 16 chirality 
patterns (xi? X 2 ) XS; X4) which are summarized in Table 

Interestingly it appears that the chirality patterns are not all independent, more precisely 
one can essentially define five classes of patterns with independent real space topology: ‘lo’i = 
{Ai, A2, Af, Ag), ^2 = (^3) ^4) A^, Ag), = (S'3, S4, 55 , Sq), = (Si, S'g), ^4^ = (S2, S7). As 
an illustration of their distinct topology, the chirality patterns associated to A3, A7, Sg, S2, Si are 
represented on Fig. [^a),..(e). Each class is characterized by the fact that it is globally invariant 
under the following three transformations: (i) chirality conjugaison C^: (xi) X2, X 3 ) X4) 
(-Xi,-X2,-X3,-X4), (ii) mirror along y-axis chains My: (xi, X2, X3, X4) (xi, X2, X4, Xs), 
(hi) translation by half unit along y-axis Ty-. (yi, X2, X35 X4) (X2, Xi) X 4 ) X 3 )- As example one 

easily verifies the identities Ag = Ty[Aj] = MylAj], Ai = CyTylAj] and A2 = C'^[A7]. Note that 
for patterns Sj in classes "Kia, '^46 one has Si = My\Si] = Ty\Si]; in that sens it is more meaningful 
to define a single class ^4 = (Si, Sg, S2, S7). In fact it then appears that the four classes are 
equivalent to classify the 16 chirality patterns in terms of the two effective chiralities 


Xa = X1X2X3X4 = sign(a2a3), yg = X3X4 = sign(bib2b3b4). 


( 4 ) 


such that % = (xa, Xb) with ^1 = (-, -), "^2 = (-, +), ^2 = (+, -) and "^4 = (+, +)• 

We terminate this section by listing other properties relating chirality patterns and LDOS. 
Pi: Since there are only 2 ‘^ patterns, it implies that there are 2 ^ (out of 2 ^) configurations of 
signs of the seven transfer integrals that leads to the same pattern (yi, y2, X 3 ) X4) and to the 
same LDOS Dj{E) (see Appendix B). In other words, without loss of generality, it is always 
possible to fix arbitrarily the sign of three out of the seven transfer integrals and hereafter the 
convention (61,62,^3) > 0 is adopted in all numerical computations. 





Table 1. Chirality patterns (xi, X2) X3) X4)- The group Ai is associated to X1X2X3X4 = “ 
whereas the group Si is associated to X1X2X3X4 = +■ 
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P2: From the properties of the moments of the LDOS (Appendix A), it is straightforward to 
show that chirality conjugation is associated to Dj{E) —)• Dj{—E). More physically it means 
that if a pattern exhibits a ZGS at 1/4 (resp. 1/2, 3 / 4 ) filling then its chirality conjugate pattern 
exhibits a ZGS at 3/4 (resp. 1/2, 1 / 4 ) filling. 

P3: In the cases ( 1 02 1 = 1 03 1 , 1 61 1 = 
pattern (xij X 2 j X35 X4)) hs My (or Ty) transformed pattern has the same DOS. 

As a consequence of properties (^2)^3), for the cases (|a2| = jasl, |6i| = I63I, I62I = |^4|) it is 
thus sufficient to study the ZGS and LDOS properties of patterns A3, A-;, Sq, S2, Si to deduce 
the properties of all other patterns. Lastly, we note that all previous numerical studies of Dirac 
points and ZGS in a-(BEDT-TTF)2l3 corresponds to the pattern A3. 


,162! = 1^41)1 using eq.(A. 6 ) it easily shown that for any 


4. Existence and stability of ZGS at 3/4 fillings 

In the following, we concentrate our numerics on the study of ZGS at 3/4 filling for patterns 
A7,A3,S'6 of classes '^i,'^2,'^3 respectively. We stress that patterns Ag,A 4,55 present very 
similar properties as A7,A3,S'6. We discard patterns of class ^4 because for the range of 
parameters that we have studied we never found ZGS at either 1/4 or 3/4 filling for these 
patterns. 


4 - 1 - ZGS and LDOS properties for homogeneous modulus |ai| = |a2| = .... = I64I = 1 
As a first illustration of the importance of the chirality pattern on the appearance of ZGS, Eig. 
[^represents respectively the LDOS Dj{E) (j = A,B,C) for patterns A3,A7,56 in the case of 
homogeneous moduli |ai| = |a2| = .... = \bi\ = 1 . On Eig. [^ for each pattern, the energies 
corresponding to 1/4 and 3/4 fillings are indicated by vertical arrows. 

Eor pattern A3, the A,B,C sites have strongly different LDOS; in addition at 1/4 and 
3/4 fillings the system is metallic, moreover there is a strong ID like VanHove singularity for 
Z>a,b,c{E) at 3/4 filling. Eor pattern A7, the LDOS verify Db{E) = Dc{E) 7^ Da{E)-, in 
addition there is a wide ZGS at 3/4 filling and a metallic state at 1/4 filling with a strong 
VanHove singularity for Da{E). Eor pattern Sq, the LDOS verify Db{E) = Dc{—E) / Da{E) 
moreover there is a narrow ZGS both for 1/4 and 3/4 fillings whereas there is a VanHove 
singularity at 1/2 filling. 

Erom these results, one can already anticipate that pattern A3 which shows a VanHove 
singulary at 3/4 filling may thus require some peculiar and probably strong anisotropy in the 
moduli in order to exhibit a ZGS state at this filling. In a similar manner, one can anticipate 
that the ZGS at 3/4 filling obtained for patterns A7 and Sq may have distinct robustness against 
the introduction of an anisotropy in the moduli. More precisely the ZGS of pattern A7 which 
exhibits a linear LDOS over wide energy range is anticipated to be more robust than the ZGS 
of pattern Sq which is relatively narrow in energy. 








Figure 2. (Color online) Real space representation of the five distinct triangular plaquette 
chirality patterns A^, ^3, Sq , Si , S'2. Triangular plaquette with Xi = + (xi = —) are in 
red (blue). Chirality conjugaison Cy- amounts to switching the color of each plaquette, in that 
way one obtains the corresponding five patterns A2, Aq, S5 , Ss , S7. By applying Ma the five 
patterns Aj, A3, Sq , Si , S2 transform into respectively ^ 48 , ^3, S5 , Si , S2- Similarly, by 
applying they transform into respectively Ag, A^, S3 , Si , S2- 


4 - 2 . ZGS and LDOS properties for BEDT moduli at pressure P = 6Kbars 

As a first illustration of the combined role of chirality pattern and inhomogeneity in 
the moduli, Fig represents the LDOS Dj{E) (j = A,B,C) for patterns A3,A7,Se 
in the case of BEDT moduli at pressure P = 6 Kbars (|a2Ma3|)|&iM^2Ufe3M^4|) 

=|ai|( 2 . 23 , 0 . 39 , 2 . 86 , 3 . 46 , 1 . 70 , 0 . 58 )[l 2 l [l 3 ], for which it is known that pattern A3 exhibits 
a ZGS at 3/4 filling. The value of |oi| is taken such that the second moment of the DOS is the 
same as for homogeneous modulus. 

At hrst sight, the main consequence of the strong inhomogeneity in the moduli is that the 
shape of the LDOS are now qualitatively much more similar for all patterns. In particular, now 
all patterns exhibit a ZGS at 3/4 filling, however for pattern A3 the ZGS is very narrow in 
energy and with strong VanHove singularities very near the Dirac point indicating a strong tilt 
of the Dirac cones. Note that for pattern Sq there is no more a ZGS at 1/4 filling where it is 
now metallic with a strong VanHove singularity. Overall, these results seem to indicate that 
chirality patterns might play a weaker role for strong inhomogeneity. 

4.3. ZGS properties for inhomogeneous moduli 

The combined roles of chirality pattern and inhomogeneity of moduli on the existence and 
stability of ZGS at 3/4 fillings are now examined more systematically for the three patterns 
A 3 ,A 7 ,S' 6 . To characterize the inhomogeneity of the moduli it is convenient to define the 


































A3 homogeneous 


A7 homogeneous 


S6 homogeneous 



Figure 3. (Color online) Local densities of states for patterns AsjAyjSe in the case of 
homogeneous modulus |ai| = \a2\ = .... = I64I = 1 , Da{E) = Da'{E) black line, Db{E) 
red line and Dc{E) blue line. 


A3 BEDT A7 BEDT S6 BEDT 



Figure 4. (Color online) Local densities of states for patterns S% in the case of BEDT 

moduli at pressure P = 6 Kbars (|a2Ma3|)|feiM^2Mfe3M^4|) =|ai|( 2 . 23 , 0 . 39 , 2 . 86 , 3 . 46 , 1 . 70 , 0 . 58 ) 
[Ill[ 13 ]. Da{E) = Da'{E) black line, Db{E) red line and Dc{E) blue line. 



























































































following six dimensionless parameters: 


_ |a2|-|ffl3l 
I«2| + |a3l’ 
_ |a2|+|a3l 

‘ 2 |ai| ’ 


A _ |fel|-|b3l 
“ |6i| + |b3p 
_ |feil + |fc3l 

“ 2 |ai| ’ 


Ab' 
Zb' ■■ 


1^21-1^41 

ifc 2 i + i^>4i ' 

1^21 + 1^41 

2 |ai| • 


(5) 


The quantities —1 < Aa^b,b' < 1 characterize the dimerization along the chains (AA’) 
and the trimerization along the two effective chains (BACA’). The ratio Za,Zb,Zb' quantify 
essentially the relative anisotropy between the different chains as compared to chain (BC). 
Within this parametrization, the case of homogeneous moduli corresponds to Aa^b,b' = 0 ^md 
Za,b,b' = 1) whereas the case of BEDT moduli at P = 6 corresponds to {za, Zb, Zb', Aq, A^, A^/) = 
(1.31,2.27,2.02,0.70,0.25,0.71). In the following, for a given set of parameters {za,b,b', Aa^fi'), 
the value of oi is determined so as to keep the second moment of the DOS equal to that 
of the homogeneous case; in that way the total bandwidth is roughly constant: |ai| = 

Y^6/[1 + z^(l + A2) + 22;^(1 + A^) + 2z^,(l + A^,)] With the convention (5i, 62 ,^ 3 ) > 0, the 
seven hopping parameters read: 


“1 = Xsloil, 

“2 = XlX3X42:a(l + Aa)|ai|, 

“3 = X2Za(l - A„)|oi|, 

bi = Zb{l + Ab)\ai\, (6) 

63 = Zb(l - Afe)|oi|, 

62 = Zb'{l + Afe/)|ai|, 

64 = XsXiZb'il - Ab/)|ai|, 

For a given chirality pattern and a given set of parameters (zafi^',Aafi^b'), to characterize the 
existence of a ZGS at 3/4 filling we proceed as follows. In a first step the energy eigenvalues 
E'n(k) are computed over a regular mesh of 400 x 400 k-points in the BZ such that the maximum 
and minimum P™“(k™“) of each band are determined. In a second step a ZGS 
at filling 3/4 is obtained whenever the following three conditions are simultaneously verfied: 

(i): / G/2; {ii): jk^^” - < J; {in): 0 < Pf“(k^“) - Pf’^’^(k§^“) < 

J <C 1 where in practice 6 ~ (vr/lOO) was taken. If either (i) or (ii) are not verified or if 

^'““(k™™) — < 0 the system is considered as metallic, whereas it is considered as 

gapped if ^““(k™™) — > 6 . Note that taking a smaller value for 6 will a priori 

decrease the number of configuration which can be considered as metallic or ZGS. 

For each of the three patterns (A3, Ay, Se), the numerical determination of possible ZGS 
at 3/4 fillings is done systematically for parameters Za^b,b' £ (0,0.2,0.4,..., 4) and for six 
distinct dimerization configurations {Aa,Ab,Ab>) as given in table such that configuration 
I is equivalent to no dimerization, II to weak and VI to strong dimerization. For a given pattern 
and a given dimerization configuration the existence of ZGS is thus explored for ( 21 )^ = 9261 
distinct configurations of parameters Za,b,b'- 

4 . 3 . 1 . histrogram distribution of the number of gapped, ZGS and metallie phases The first 
statistic FigJ^ represents the histogram distribution of the number of gapped, ZGS and metallic 
configurations for patterns {A 3 ,Ay,SQ) at 3/4 filling and for the six distinct dimerization 
configurations (I,II,...VI). 

The common and expected feature is that the number of gapped configurations (black bars) 
increases with the strength of the dimerization along the different chains; for all patterns it 
eventually becomes the dominant phase for strong dimerization (VI). Despite this common 
feature, we find it relatively surprising that for pattern Sq, the number of gapped configurations 
is always relatively large even for low or no dimerization. As far as the number of metallic 









Table 2. Dimerization configurations (A^, Af,, A;,/) 



I 

II 

III 

IV 

V 

VI 

Aa 

0 

0.25 

0.25 

0.5 

0.75 

0.75 

Ab 

0 

0.25 

0.75 

0.5 

0.25 

0.75 

Ay 

0 

0.25 

0.25 

0.5 

0.75 

0.75 


configurations (blue bars) is concerned, its behaviour is very different for the three distinct 
patterns {A 3 , Ay, Sq). For pattern ^3, it represents the dominant phase for small dimerization 
(I) and then it steadily decreases with increasing dimerization and it becomes the rarest phase for 
large dimerization (VI); this type behaviour is somehow as expected. For pattern Aj, whatever 
the dimerization, the number of metallic configurations is surprisingly always very small in 
proportion such that roughly speaking this pattern exhibits either a ZGS or a gapped phase. 
For pattern Sq, the number of metallic configurations slowly decreases when dimerization is 
increased and it is always below that of gapped and ZGS phases excepted for dimerization 
(III) where it represents the dominant phase. Overall, we find it quite remarkable that for 
most dimerization configurations and for the three patterns (^3, Ay,Sq), the number of metallic 
configurations remains always smaller than the number of ZGS phases. In fact, regarding the 
number of ZGS configurations, the three patterns {A 3 , Aj, Sq) show distinctive behaviours. For 
patterns {A 3 ,Sq) the number of ZGS configurations shows no clear tendency as a function of 
dimerization strength; for each dimerization configuration it constitutes around 1/3 of the 9261 
anisotropy configurations, slightly higher for Sq and slightly lower for A 3 . More remarkably, for 
pattern Ay, even if the number of ZGS configurations decreases steadily when the dimerization 
strength is increased, it still constitutes the largely dominant phase over the entire studied range 
of dimerization. In summary, the main message to retain from Figj^ is that over the range of 
anisotropy and dimerization that we have considered, the pattern Ay is by far the most favorable 
for the appearance of stable ZGS at 3/4 filling whereas patterns A 3 , Sq are more or less similarly 
favorable. 

4-3.2. Anisotropy versus dimerization dependency of ZGS properties The previous histograms 
mainly provide quantitative information of the role played by the dimerization parameters 
(Aa, Afc, Afc/) on the appearance of ZGS at 3/4 filling, however they do not allow to understand 
the role played by the anisotropy parameters Za^b,b'- In order to quantify the role played by the 
anisotropy between the different chains, for each dimerization configuration we have computed 
both the average value and the variance of the anisotropy for configurations Za b,b' that lead to 
ZGS at 3/4 filling for patterns A 3 , Ay, Sq. The results are summarized on Fig. On this figure 
the grey region corresponds to the case of a uniform distribution of ZGS configurations over the 
explored interval of anisotropy : {za, Zb, zy) G [0,4]; we reminds that no anisotropy correspond to 
Za = Zb = zy = ^. For pattern A 3 , for no or weak dimerization there is a strong deviation from 
the uniform distribution, more precisely the more favorable anisotropy configurations correspond 
to large value of Za and small value of Zb, zy. By increasing the dimerization the average values of 
Za, Zb, zy enter the region corresponding to a uniform distribution. For pattern Sq the behaviour 
is somewhat symmetric to that of pattern A 3 , for small dimerization small value of Za and larger 
value of Zb,zy are more favorable to obtain a ZGS at 3/4 filling. Finally for pattern Ay the 
distribution of ZGS is almost independent on both the anisotropy and dimerization. The main 
message to retain from this Fig. is that the larger the dimerization the less important is the 
role of the anisotropy; in other words the strength of the dimerization along the different chains 
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A7, 3/4 filling 
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Figure 5. histogram distribution of the number of configurations that are gapped (black), 
ZGS (red) or metallic (blue) for patterns ^ 3 , A^, Sq at 3/4 filling for the six distinct dimerization 
configurations (I,II,...VI). 


seems a more pertinent parameter than the anisotropy between the chains. 

4 . 3 . 3 . k-space distribution of Dirac points associated to ZGS phases As a last characterization 
of the ZGS properties, the following Fig. represents the evolution of the k-space distribution 
of the Dirac points associated to the ZGS at 3/4 for pattern A 3 as function of the dimerization. 
Figs. |8|9| represent this same quantity but for pattern A 7 and Sq respectively. The most striking 
feature on these figures is that the case of no dimerization (I) is very different from all the other 
cases. In fact, in the absence of dimerization the distribution of Dirac points associated to ZGS 
is much less spreaded over the Brillouin zone, moreover it differs strongly from one pattern to the 
other. By increasing the dimerization, the k-space distribution of ZGS’s Dirac points for patterns 
A 3 and Ay become more and more similar, in fact they seem to only differ by number of Dirac 
points (intensity) as already seen in their corresponding histogram distribution. For sufficiently 
large dimerization (V,VI), the k-space distribution of ZGS’s Dirac points for all patterns are 
very similar. All these features are further indications that the larger the dimerization the less 
important is the role played by the form of the chirality pattern. 
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A3, 3/4 filling 


A7, 3/4 filling 
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Figure 6. (color online) Averaged anisotropy parameter values Za (red), zi, (black) zi,/ (blue) 
versus dimerization configurations (I,II,...VI) for the ZGS configurations found for patterns A 3 , 
Aj, S 3 at 3/4 filling. The errorbars encode the meansquare fluctuations. For comparison, the 
grey region corresponds to the case of a uniform distribution of ZGS configurations over the 
interval {za,Zb,Zb') G [0,4]. 



Figure 7. (color online) k-space distribution of the Dirac points associated to the ZGS 
configurations obtained for patterns A 3 at 3/4 filling for the six distinct dimerization 
configurations (I,II,...VI). 











































Figure 8 . k-space distribution of the Dirac points associated to the ZGS configurations 
obtained for patterns at 3/4 filling for the six distinct dimerization configurations (I,II,...VI). 



Figure 9. k-space distribution of the Dirac points associated to the ZGS configurations 
obtained for patterns S% at 3/4 filling for the six distinct dimerization configurations (I,II,...VI). 














5. Summary 

In summary, we have examined the existence of ZGS at 3/4 filling for a tight binding model with 
seven nearest neighbors transfer integrals (ai, 02 , 03 , 6i, 62 ; ^ 4 ) between the four molecules of 

the unit cell, as in a-type organic conductor. We have systematically studied the combined role 
of the sign and moduli of the transfer integrals on the appearance of ZGS at 3/4 filling. We have 
shown that for given moduli, the are only 16 independent band structures Ai, ■ ■ ■ ,Ss. Each band 
structure is characterized by a chirality pattern {xi,X2-, XS; X 4 )) where each chirality encodes the 
sign of the product of three transfer integrals that compose one of the four distinct triangular 
plaquette of the unit cell (xi = sign(a26i64), X2 = sign(a36263), X3 = sign(ai6i62), X4 = 
sign(ai 6364 )). In term of this classification, it appears that all previous works on q;-(BEDT- 
TTF) 2 l 3 corresponds to pattern A 3 with (xi, X2, X3, X4) = 

As far as the properties of the LDOS and ZGS at 3/4 (resp. 1/4) filling are concerned, 
the 16 patterns can be further classified into only four different classes determined by two 
effective chiralities Xa = sign(aia 2 ) and Xb = sign( 6 i 626364 ). Using extensive numerics we have 
established that for both large, or no inhomogeneity in the moduli |ai|, ...I 64 I, the patterns A 7 , Ag 
(resp. Ai, A 2 ) of class (xa, Xb) = are the most favorable patterns for the existence of ZGS 

at 3/4 (resp. 1/4) filling. At the opposite, for patterns {Si, S2, Sj, Sg) of class {Xa,Xb) =(+;+) 
we never found any ZGS at either 1/4 or 3/4 filling. Finally, patterns (S' 5 , Sq) (resp. S3, S' 4 ) and 
A 3 ,Aq (resp. ^ 3 ,^ 4 ) of class {Xa,Xb) = (+> ~) and (Xa,Xb) = (~)+) present an intermediate 
situation; they allow for ZGS at 3/4 (resp. 1/4) filling but are much less favorable that patterns 
of class {Xa,Xb) = (—)—)• We note however that the role played by the chirality pattern is 
dependent on strength of the dimerization in the different effective chains formed by the transfer 
integrals. More precisely for weak or no dimerization the form of the chirality pattern strongly 
determines the probability to find a ZGS at 1/4 or 3/4 filling whereas for strong dimerization the 
possibility of having a ZGS become more or less equivalent for the three classes (—,—),(—, +) 
and (+, —). From that perspective we note that the usual transfer integral parameters taken to 
modelize Q;-(BEDT-TTF) 2 l 3 are in the strong dimerization regime. 

On a larger perspective, we note that for bipartite lattice (like square lattice) when all 
plaquette chiralities are negative it corresponds to the so called uniform vr-flux phase for which it 
is well known that there is a ZGS at 1/2 filling. Our results suggest that on a non-bipartite lattice 
like the triangular lattice, it is necessary to have a fractional density p/q of vr-flux plaquettes to 
obtain a ZGS at some filling r/g; the difficulty in such a situation is that there appears many 
non equivalent chirality patterns. 
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Appendix A. Local density of states and chirality patterns 

Consider the a-BEDT tight binding model with only nearest-neighbor hopping on the triangular 

I O •j-'J-'l 

^ = EE[%i*)oi+h.c.], (A.i) 

i <j,i> 

with tij G (ai,..., 64 ). The one particle local Greens function on a site j is Gjj{E) = 
{j\{E — H)~^\j). It can be expanded as 




1 {j\H^\j) 

E ^ E^ 

t 


(A.2) 



The local density of states is obtained from the identity Dj{E) = + iO)}, which 

allows to write 

/ oo 

dE E^D,{E) = {j\H%)= Y, E - E (A.3) 


This last identity means that the moment of the local density of states on site j is equal 
to sum over all closed path of steps that start from j and return to j. The key property is that 
each path contributes with a weight proportionnal to the product of moduli of all the hoppings 
along the path and with a sign given by the product of the chiralities of all the triangular 
plaquettes encircled by the path. For the tight binding model of a-BEDT, the moment for each 
site up to ^ = 3 are easily obtained as: 


T 03 + ^3 + ^4, 

= 2(04 + 62 + ^ 3 )) 

= 2 {af + + 64 ), 

= 2(y2|a3fe2^3| + Xl|o2&1^4|) + X3|aifel^2| + X4|oi&3^4|, 
= 2(X2|a362ft3| + X3|aiftl^2| + X4|ai63^4|), 

Ji}, = 2(xi|a26i64| + X3|ai^i&2| + X4|ai&3&4|)- 

More generally, odd and even moment have the generic form: 

= EaXa XaX/3X7 4/37,i> 

= V] + XaX/3 + XlX2X3X4 


(A.4) 


(A.5) 


where > 0 depend only on the moduli (|ai|,| 64 |)- For 1 02 ! = |a3|, 

1^1 1 = 1^3! and 1621 = 1^41, one easily obtains that the global moments 

verify 

^2£+i ^ X2)(x^ + X3X4a:f34) + (X3 + X4)(a:3 + XiX2a:i23), gN 

+ (xi+ X2)(x3+ X4)yi3+ xiX 2 yi 2 + X3X4y34+ X1X2X3X4 yW- ^ 


Appendix B. The 8 transfer integral configuration associated to pattern A 3 


Table Bl. The 8 = 2^ configurations of signs of the seven transfer integrals that lead to the 
chirality pattern A^. All these configurations have the same LDOS Dj[E). Strictly speaking 
for k-space dependent properties it is still necessary to distinguish two groups, i.e configurations 
(1,2,5, 6 ) with (6164 < 0 and 6263 < 0) and congigurations (3,4,7, 8 ) with (6164 > 0 and 62^3 > 0). 
The energy bands -En(k) of the first group become equal to that of the second group after the 
substitution kx ^ kx + vr. 


^3 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

sign(ai) 

+ 

+ 

+ 

+ 

- 

- 

- 

- 

sign(a2) 

+ 

+ 

- 

- 

+ 

+ 

- 

- 

sign(a3) 

- 

- 

+ 

+ 

- 

- 

+ 

+ 

sign(6i) 

+ 

- 

+ 

- 

+ 

- 

+ 

- 

sign( 62 ) 

- 

+ 

- 

+ 

+ 
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+ 
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sign( 63 ) 

+ 

- 

- 

+ 

- 

+ 

+ 

- 

sign( 64 ) 

- 

+ 

+ 

- 

- 

+ 

+ 
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